Drawing with Complex Numbers 



Michael Eastwood^ and Roger Penrose 

It is not commonly realized that the algebra of complex numbers can be 
used in an elegant way to represent the images of ordinary 3-dimensional 
figures, orthographically projected to the plane. We describe these ideas 
here, both using simple geometry and setting them in a broader context. 

Consider orthogonal projection in Euclidean n-space onto an m-dimen- 
sional subspace. We may as well choose coordinates so that this is the 
standard projection P : W 1 — > IR m onto the first m variables. Fix a non- 
degenerate simplex E in M". Two such simplices are said to be similar if 
one can be obtained from the other by a Euclidean motion together with 
an overall scaling. This article answers the following question. Given n + 1 
points in M m , when can these points be obtained as the images under P of 
the vertices of a simplex similar to E? 

When n = 3 and m = 2, then P is the standard orthographic projection 
(as often used in engineering drawing) and we are concerned with how to 
draw a given tetrahedron. We shall show, for example, that four points 
a,/3,7, 8 in the plane are the orthographic projections of the vertices of a 
regular tetrahedron if and only if 

(a + (3 + 1 + 5) 2 = A(a 2 + (3 2 + 1 2 + 5 2 ) (1) 

where a, (3, 7, 5 are regarded as complex numbers! Similarly, suppose a cube 
is orthographically projected and normalised so that a particular vertex is 
mapped to the origin. If a, (3, 7 are the images of the three neighbouring 
vertices, then 

« 2 + / 3 2 + 7 2 = 0, (2) 

again as a complex equation. Conversely, if this equation is satisfied, then 
one can find a cube whose orthographic image is given in this way. Since 
parallel lines are seen as parallel in the drawing, equation allows one to 
draw the general cube: 
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In this example, a = 2 — 2Gi 
(3 = -23 + 2i 
7 = 14 + 7i 








The result for a cube is known as Gauss' fundamental theorem of 
axonometry — see [§, p. 309] where it is stated without proof. In engineer- 
ing drawing, one usually fixes three principal axes in Euclidean three-space 
and then an orthographic projection onto a plane transverse to these axes 
is known as an axonometric projection (see, for example, || Chapter 17]). 
Gauss' theorem may be regarded as determining the degree of foreshortening 
along the principal axes for a general axonometric projection. The projec- 
tion corresponding to taking a, (3, 7 to be the three cube roots of unity is 
called isometric projection because the foreshortening is the same for the 
three principal axes. In an axonometric drawing, it is conventional to take 
the image axes at mutually obtuse angles: 



If I a I = a, \/3\ — b, I7I = c, then equation (§) is equivalent to the sine rule 
for the triangle with sides a 2 , (3 2 , j 2 , namely 



In this form, the fundamental theorem of axonometry is due to Weisbach and 
was published in Tubingen in 1844 in the Polytechnische Mitteilungen of Volz 
and Karmasch. Equivalent statements can be found in modern engineering 
drawing texts (e.g. 0, p. 44]). 
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Equation (g) may be used to give a ruler and compass construction of 
the general orthographic image of a cube. If we suppose that the image of a 
vertex and two of its neighbours are already specified, then (Q) determines (up 
to a two-fold ambiguity) the image of the third neighbour. The construction 
is straightforward except perhaps for the construction of a complex square 
root for which we advocate the following as quite efficient: 




Firstly ( is constructed by marking the real axis at a distance \\z\\ from the 
origin. Then, a circle is constructed passing through the three points (, 1, 
and z. Finally, the angle between 1 and z is bisected and yfz appears where 
this bisector meets the circle. 

In engineering drawing, it is more usual that the images of the three 
principal axes are prescribed or chosen by the designer and one needs to 
determine the relative degree of foreshortening along these axes. There is a 
ruler and compass construction given by T. Schmid in 1922 (see, for example, 
[|, §17.17-17.19]): 
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In this diagram, the three principal axes and a are given. By drawing a 
perpendicular from a to one of the of the principal axes and marking its 
intersection with the remaining principal axis, we obtain P. The point Q 
is obtained by drawing a semi-circle as illustrated. The point R is on the 
resulting line and equidistant with a from Q. Finally, 3 is obtained by 
dropping a perpendicular as shown. It is easy to see that this construction has 
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the desired effect — in Euclidean three-space rotate the right-angled triangle 
with hypotenuse Pa about this hypotenuse until the point Q lies directly 
above in which case R will lie directly above (3 and the third vertex will lie 
somewhere over the line through and Q. One may verify the appropriate 
part of Weisbach's condition 



sin2v4 sin 25 



(3) 



by the following calculation. Without loss of generality we may represent 
all these points by complex numbers normalised so that Q — 1. Then it is 
straightforward to check that 

„ „ a(a + a) + 2(1 — a — a) „ a(a + a) + 2(1 — a — a) . 
R=l + t-ta, P= — - ^ -, (3=— 4- -% 



a — a 



2 — a — a 



and therefore that 



a 2 + (3 2 



[a 



l)(a - l)(a + a - 1) 



(a + a- 2) : 



That a 2 + [3 2 is real is equivalent to @. 

To prove Gauss' theorem more directly consider three vectors in M 3 as the 
columns of a 3 x 3 matrix. This matrix is orthogonal if and only if the three 
vectors are orthonormal. It is equivalent to demand that the three rows be 
orthonormal. However, any two orthonormal vectors in M 3 may be extended 
to an orthonormal basis. Thus, the condition that three vectors 



X\ 










k 2/1 > 




. 2/2 , 
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in M 2 be the images under P : M. 3 — > IR 2 of an orthonormal basis of M. 3 , is 
that 

( xi x 2 x 3 j and ^ y x y 2 y 3 j 

be orthonormal in R 3 . Dropping the overall scale, we obtain 

x 2 + x 2 2 + x 3 2 = y 2 + y 2 + y 3 2 and x x y x + x 2 y 2 + x 3 y 3 = 0. 

Writing, a = Xi+iyi, [3 = x 2 +y 2 , 7 = x 3 +y 3 , these two equations are the real 
and imaginary parts of (0). To deduce the case of a regular tetrahedron as 
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described by equation ([I]) from the case of a cube as described by equation (|2|), 
it suffices to note that equation ([!]) is translation invariant and that a regular 
tetrahedron may be inscribed in a cube. Thus, we may take 5 = a + [3 + 7 
and observe that (|]) and (|2|) are then equivalent. 

It is easy to see that the possible images of a particular tetrahedron X in 
M. 3 under an arbitrary Euclidean motion followed by the projection P form 
a 5-dimensional space — the group of Euclidean motions is 6-dimensional but 
translation orthogonal to the plane leaves the image unaltered. It therefore 
has codimension 3 in the 8-dimensional space of all tetrahedral images (2 de- 
grees of freedom for each vertex). Allowing similar tetrahedra rather than 
congruent reduces the codimension to 2. Therefore, two real equations are to 
be expected. Always, these two real equations combine as a single complex 
equation such as ([!]) or (0). At first sight, this is perhaps surprising and even 
more so when the same phenomenon occurs for P : M. n — > M 2 for arbitrary n. 

For n = 3, there is a proof of Gauss' theorem which brings in complex 
numbers at the outset. Consider the space H of Hermitian 2x2 matrices 
with zero trace, i.e. matrices of the form 



X 



w 

u — iv 



u + IV 





' u ' 


for 


V 




k w > 



We may identify H with R 3 and, in so doing, — det X becomes the square of 
the Euclidean length. The group G of invertible 2x2 complex matrices of 
the form 

[b a J 



A 



acts linearly on H by X 1— > AXA 

det(AXA*) : 



Moreover, 
= (|a| 2 + |6| 2 ) 2 detX 



so G acts by similarities. It is easy to check that all similarities may be 
obtained in this way. (This trick is essentially as used in Hamilton's theory 
of quaternions and is well known to physicists — in modern parlance it is 
equivalent to the isomorphism of Lie groups Spin(3) = SU(2).) Therefore, 
an arbitrary orthographicimage of a cube may be obtained by acting with A 
on the standard basis 
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and then picking out the top right hand entries. We obtain 



A 



A 



A 
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* i(a 2 + b 2 ) 

* * 
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i(a 2 + b 2 
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and therefore a 2 + (3 2 + •y 2 = 0, as required. Conversely, this is exactly the 
condition that a, (3, 7 may be written in this form. (Compare the half angle 
formulae— if s 2 + c 2 = 1, then s = 2t/(l + t 2 ) and c = (1 - t 2 )/(l + t 2 ) for 
some t.) That Gauss [|3], p. 309] makes the same observation concerning the 
form of a, [3, 7 suggests that perhaps he also had this reasoning in mind. 

The proof of Gauss' theorem using orthogonal matrices clearly extends 
to P : M n — > 1R 2 = C for arbitrary n. To state it, the following terminology 
concerning the standard projection P : W 1 — > M m is useful. We shall say 
that t>i,t>2, ■ ■ ■ ,v n G M m are normalised eutactic if and only if there is an 
orthonormal basis u\, u%, . . . , u n of M. n with Vj = Puj for j = 1,2, ... ,n. We 

are eutactic if and only if fivi, \iVi, . . . , fiv n 



shall say that V\, v%, ■ ■ ■ , v n G 
are normalised eutactic for some \x 7^ 0. 



Theorem TTie points z\,Z2, 



Zn G 



are eutactic if and only if 



zi 2 + z 2 2 + 



and not all Zj are zero. 

There is a proof for n = 4 based on the isomorphism 

Spin(4) SU(2) x SU(2) 

and, indeed, this is how we came across the theorem in the first place. How- 
ever, a more direct route to complex numbers and one which applies in all 
dimensions is based on the observation that Gr^~(IR 2 ), the Grassmannian of 
oriented two-planes in M n , is naturally a complex manifold. When n — 3, 
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this Grassmannian is just the two-sphere and has a complex structure as the 
Riemann sphere. In general, consider the mapping 

CP n _i \ RP n _! Gr+(R n ) 

induced by C n 3 z i— » izAz. In other words, a complex vector z = x+iy G C™ 
is mapped to the two-dimensional oriented subspace of lR n spanned by x and 
y, the real and imaginary parts of z. Let ( , ) denote the standard inner 
product on M n extended to C n as a complex bilinear form. Then, (z, z) = 
imposes two real equations 

ll x ll 2 = \\y\\ 2 an d ( x ,y) = Q 

on the real and imaginary parts. In other words, x, y is proportional to an or- 
thonormal basis for span{x, y}. Hence, if z and w satisfy (z, z) = = (w, w) 
and define the same oriented two-plane, then w = Xz for some A G C \ {0}. 
The non-singular complex quadric 

K = {[z]e CP n _i s.t. (z,z) = 0} 

avoids PJP n _i C CP n _i and we have shown that h\k is injective. It is clearly 
surjective. The isomorphism 

Ti : K — — > Gr+(R n ) 

respects the natural action of SO(n) on K and Gr^(IR n ). The generalised 
Gauss theorem follows immediately since, rather than asking about the image 
of a general orthonormal basis under the standard projection P : W 1 — > M 2 , 
we may, equivalently, ask about the image of the standard basis ei, e 2 , . . . , e n 
under a general orthogonal projection onto an oriented two-plane II C M n . 
Any such II is naturally complex, the action of i being given by rotation by 
90° in the positive sense. If II is represented by [zi, Z2, ■ ■ ■ , z n ) £ K as above 
and we use x, y G II to identify II with C, then ej t— > Zj and 

Zl 2 + z 2 2 + ■ ■ ■ + z n 2 = (z, z) = 0, 

as required. Conversely, a solution of this complex equation determines an 
appropriate plane IT. 

For the case of a general tetrahedron or simplex and for general m and 
n, it is more convenient to start with Hadwiger's theorem Q or [0, page 251] 
as follows. The proof is obtained by extending our orthogonal matrix proof 
of Gauss' theorem. 
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Theorem (Hadwiger) Assemble i>i,i>2, ,v n G M m as the columns of an 
mxn matrix V. These vectors are normalised eutactic if and only ifVV 1 = 1 
(the m x m identity matrix). 

Proof If v i, i>2, • • • , v n are normalised eutactic, then assembling a corre- 
sponding orthonormal basis of M ra as the columns of an n x n matrix, we 
have V = PU and U t U = 1 (the nxn identity matrix). Therefore, UU l = 1 
and 

vv l = PUU 1 ? 1 = PP l = 1, 

as required. Conversely, if VV* = 1, then the columns of V 1 may be com- 
pleted to an orthonormal basis of R n , i.e. V 1 = U t P t for UU* = 1. Now, 
U l U = 1 and V = PU, as required. □ 

The case of a general simplex is obtained essentially by a change of basis 
as follows. Suppose ai, a 2 , . . . , a n , a n+ i are the vertices of a non-degenerate 
simplex £ in IR n whose centre of mass is at the origin. In other words, the 
n x (n + 1) matrix A has rank n and Ae = where e is the column vector all 
of whose n + 1 entries are 1. Form the (n + 1) x (n + 1) symmetric matrix 

Q = A*(A4*)- 2 A, 

noting that rank A = n implies the moment matrix AA l is invertible. 

Theorem Given bi, 62, • • • , b n , b n+ \ G M. m assembled as the columns of an 
mx (n+1) matrix B, these vectors are the images under orthogonal projection 
P : W l — > R m 0/ £/ie vertices of a simplex congruent to S if and only if 

BQ l B = 1. (4) 

Proof The vertices of a simplex congruent to E are the columns of a matrix 
UA + ae l for some orthogonal matrix U and translation vector a G M™. Also, 
note that Qe = 0. Thus, if B = P(UA + ae*), then 

BQB 1 = PUAQAWP 1 

= PU AA l (AA t )~ 2 AA t U t P t 
= PUU l P l = PP l = 1, 
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as required. Conversely, Qe = implies that (^) is translation invariant. So, 
without loss of generality, we may suppose that b\ + b 2 + ■ ■ ■ + b n + b n+ i = 0, 
that is to say, Be = 0. Writing out (|) in full gives 

BA t (AA t y 1 (BA t (AA t )~ 1 ) t = 1 

so, by Hadwiger's theorem, there is an orthogonal matrix U so that 

BA t (AA t y 1 = PU. 

Thus, 

BA t (AA t )~ 1 A = PUA and Be = 0. 

Certainly, B = PUA is a solution of these equations but it is the only 
solution since A t (AA t )~~ 1 A has rank n and e is not in the range of this linear 
transformation. □ 



Corollary (case m = 2) Points Zi, z 2 , . . . , z n , z n+ i G C are the images un- 
der orthogonal projection of the vertices of a simplex similar to E if and only 
if 

z f Qz = 

where z is the column vector with components z±, z 2 , ■ ■ ■ , z n , z n+ \. 

It is, of course, possible to compute Q explicitly for any given example. If 
the simplex S has some degree of symmetry, however, we can often circum- 
vent such computation. Consider, for example, the case of a regular simplex. 
From the corollary above, we know that the image of such a simplex in the 
plane is characterised by a complex homogeneous quadratic polynomial. The 
symmetries of the regular simplex ensure that this polynomial must be in- 
variant under S n+ i, the symmetric group onn + 1 letters. Hence, it must be 
expressible in terms of the elementary symmetric polynomials. Equivalently, 
it must be a linear combination of 

(zi + z 2 H h z n + z n+ i) 2 and z\ + z 2 2 H h z n 2 + z n+ i 2 . 

Up to scale, there is only one such combination which is translation invariant, 
namely 

(z 1 + z 2 H h z n + z n+1 ) 2 -(n+ l)(zi 2 + z 2 H h z 2 + z n+1 2 ). (5) 
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It follows that the vanishing of this polynomial is an equation which charac- 
terises the possible images of a regular simplex under orthogonal projection 
into the plane. The special case n = 2 characterises the equilateral triangles 
in the plane [p], Problem 15 on page 79]. 

Equation (0) characterising the orthographic images of a cube, may be 
deduced by similar symmetry considerations. If a particular vertex is mapped 
to the origin and its neighbours are mapped to a, (3, 7 then, since each of these 
neighbouring vertices is on an equal footing, the polynomial in question must 
be a linear combination of (a + (3 + 7) 2 and a 2 + f3 2 + 7 2 . To find out which 
linear combination we need only consider a particular projection, for example: 

7 = i 



a = (3 = 1 

In this example, (a + (3 + •y) 2 = 2i and a 2 + /3 2 + 7 2 = 0. Up to scale, 
therefore, (Q) is the correct equation. 

The case of a regular dodecahedron is similar. Using the fact that a cube 
may be inscribed in such a dodecahedron J5|, we may deduce a particular 
projection: 




with (a + 13 + 7 ) 2 = (7 - 3^5)/2 and a 2 + (3 2 + 7 2 = (2 - v / 5)/2. In this 
particular case, 

(a + (3 + 7 ) 2 + (VE - l)(a 2 + (3 2 + 7 2 ) = 0. 
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Therefore, this is the correct equation in the general case. It may be used 
as the basis of a ruler and compass construction of the general orthographic 
projection of a regular dodecahedron. 

It is interesting to note that if all the vertices of a Platonic solid are 
orthographically projected to z%, z 2 , ■ ■ ■ , Zn £ C, then 

(zi + z 2 + • ■ • + z N f = N( Zl 2 + z 2 2 + ■■■ + z N 2 ) (6) 

(compare ([5])). For a tetrahedron, this is just equation ([[]). To verify 
(||) for the other Platonic solids, firstly note that it is translation invari- 
ant. Therefore, it suffices to impose z\ + z 2 + ■ • ■ + £jv = and show that 
Zi 2 + z 2 2 + - ■ - + z N 2 = 0. The case of a cube now follows immediately since its 
vertices may be grouped as two regular tetrahedra. The dodecahedral case 
may be dealt with by grouping its vertices into five regular tetrahedra. The 
regular octahedron is amenable to a similar trick but not the icosahedron. 
Rather than resorting to direct computation, a uniform proof may be given as 
follows. As before, assemble the vertices of the given solid £ as the columns 
of a matrix A, now of size 3 x N, and consider the moment matrix M = AA t . 
Observe that 



(no) 



M 



1 

i 




Zl 2 + Z 2 2 + ■ ■ ■ + 27V 2 - 



The moment matrix is positive definite and symmetric. In other words, it 
defines a metric on M 3 , manifestly invariant under the symmetries of S. If S 
is regular — or, more generally, enjoys the symmetries of a regular solid (e.g. a 
cuboctahedron or rhombicosidodecahedron) — then its symmetry group acts 
irreducibly on M. 3 . Thus, M must be proportional to the identity matrix and 
the result follows. For a general solid S, the two complex numbers 

±vV + 2 2 2 + ... + 2 "2 

are the foci of the ellipse 

where R is the inverse of the quadratic form obtained by restricting M to the 
plane of projection. 
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This reasoning also works in higher dimensions where it shows (as con- 
jectured to us by H.S.M. Coxeter) that the orthogonally projected images 
in the plane of the N vertices of any regular polytope, real or complex, will 
satisfy equation @. Of course, this excludes regular polygons (whose sym- 
metry groups act reducibly except in dimension two) orthographic images of 
which will satisfy @ if and only if the image is itself regular. For polyhedra 
other than simplices, a quadratic equation such as (||) is no longer sufficient 
to characterise the orthogonal image up to scale. In general, there will also 
be some linear relations. For a non-degenerate iV-gon there will be N — n — 1 
such relations. The simplest example is a square in M? which is characterised 
by the complex equations 

(a + (3 + 7 + 5f = 4(a 2 + (3 2 + 7 2 + 5 2 ) and a + 7 = f3 + 5. 

It is interesting to investigate further the relationship between a non- 
degenerate simplex £ in IR n and its quadratic form Q = A t (AA t )~ 2 A. Recall 
that A is the n x (n + 1) matrix whose columns are the vertices of S. There 
are several other formulae for or characterisations of Q. Let S denote the 
(n + 1) x (n + 1) symmetric matrix 



1- 



n+1 



1 1 
1 1 



1 1 



It is the matrix of orthogonal projection in IR n+1 in the direction of the 
vector e. We maintain that Q is characterised by the equations 

QA*A = S and Qe = 0. 

Certainly, if these equations hold, then they are enough to determine Q 
because the moment matrix M = A 1 A has rank n and e is not in its range. 
The second equation is evident and the first equation with Q replaced by 
A t (AA t )~ 2 A and simplified reads 

A t (AA t y 1 A = S. 

To see that this holds it suffices to observe that it is clearly true after post- 
multiplication by A t or e. We may equally well characterise Q by means of 
the equations 

A f AQ = S and Qe = 0. 
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These equations relate M and Q geometrically — both matrices annihilate e 
whilst on the hyperplane orthogonal to e they are mutually inverse. This 
implies that M and Q are generalised inverses || of each other. Thus, 

Q M {A 1 A) j = A ] A ]t 

where A* is the generalised inverse of A. In this case, A* = A t (AA t )~ 1 . This 
also shows how to compute Q more directly in certain cases. The moment 
matrix M has direct geometric interpretation as the various inner products of 
the vectors ai, a 2 , . . . , a n , a n+ \. In the case of a regular simplex, for example, 
we know that ||aj|| 2 is independent of i, that ||aj — aj\\ 2 is independent of 
i ^ j, and that a\ + a 2 + . . . + a n + a n+ i = 0. We may deduce that, with a 
suitable overall scale, M — S. Since S' = S, it follows that Q = S. This is 
a direct derivation of (||). 

It is clear geometrically that M or, equivalently, Q determines S up to 
congruency. Alternatively, one can argue algebraically — it is easy to check 
that if A 1 A = B l B, then U = AA^BA 1 )' 1 is orthogonal and A = UB. 
Therefore, the possible quadratic forms Q which can arise give a natural 
parametrisation of the non-degenerate simplices up to congruency. Choosing 
a basepoint So with corresponding matrix Aq, and mapping X e GL(n, M) 
to X _1 y4 identifies the space of non-degenerate simplices up to congruency 
with the homogeneous space GL(n, R)/0(n). This homogeneous space may 
also be identified with the space of positive definite n x n quadratic forms 
by sending X G GL(n, M) to XX 1 . The (n+ 1) x (n+ 1) quadratic form Q 
corresponding to X~ 1 Ao is given by AqXX 1 Aq . It follows that the general 
Q which can arise is characterised by the following two conditions: 

• Qe = and only multiples of e are in the kernel of Q. 

• All other eigenvalues of Q are positive. 

It is also possible to repeat this analysis in pseudo-Euclidean spaces. The 
only difference is that the condition that the non-zero eigenvalues of Q be 
positive is replaced by a condition on sign precisely reflecting the original 
signature of the inner product. 

Finally we should mention some possible applications. There is much 
current interest in computer vision. In particular, there is the problem of 
recognising a wire-frame object from its orthographic image. The results we 
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have described can be used as test on such an image, for example to see 
whether a given image could be that of a cube or to keep track of a moving 
shape. It is clear that such tests could be implemented quite efficiently. An- 
other possibility is in the manipulation of CADDfl data. Rather than storing 
an image as an array of vectors in M 3 , it may be sometimes be more efficient to 
store certain tetrahedra within such an image by means of the corresponding 
quadratic form. For orthographic imaging this may be preferable. 

We would like to thank H.S.M. Coxeter for drawing our attention to 
Hadwiger's article, R. Michaels and J. Cofman for pointing out Gauss' and 
Weisbach's work, and E.J. Pitman for many useful conversations. 
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